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The recent discovery of ten new dwarf galaxy candidates by the Dark Energy Survey (DES) and the 
Panoramic Survey Telescope and Rapid Response System (Pan-STARRS) could increase the Eermi 
Gamma-Ray Space Telescope’s sensitivity to annihilating dark matter particles, potentially enabling 
a definitive test of the dark matter interpretation of the long-standing Galactic Genter gamma-ray 
excess. In this paper, we compare the previous analyses of Eermi data from the directions of the 
new dwarf candidates (including the relatively nearby Reticulum II) and perform our own analysis, 
with the goal of establishing the statistical significance of any gamma-ray signal from these sources. 

We confirm the presence of an excess from Reticulum II, with a spectral shape that is compatible 
with the Galactic Genter signal. The significance of this emission is greater than that observed from 
99.84% of randomly chosen high-latitude blank-sky locations, corresponding to a local detection 
significance of 3.2cr. We improve upon the standard blank-sky calibration approach through the use 
of multi-wavelength catalogs, which allow us to avoid regions that are likely to contain unresolved 
gamma-ray sources. 


I. INTRODUCTION 

Over the past several years, a bright and statisti¬ 
cally significant excess of gamma-rays has been reported 
from the region surrounding the Galactic Center 
[9]. The spectral and morphological characteristics of 
this signal are each in good agreement with that pre¬ 
dicted from annihilating dark matter particles with a 
mass of m-oM ^ 35-60 GeV and a cross section of av ^ 
10“^^ cm^s“^ (for the representative case of annihilations 
to bb). And although the proposed astrophysical explana¬ 
tions for this excess have been shown to face considerable 
challenges, it is not currently possible to entirely rule out 
the possibility that these photons originate from a large 
population of unresolved point sources nnHH] or from a 
series of cosmic ray outbursts [HHIS]- In light of this situ¬ 
ation, gamma-ray observations of the Milky Way’s dwarf 
spheroidal galaxies play a critical role, being potentially 
able to provide a confirmation or refutation of the dark 
matter interpretation of the Galactic Center excess. 

Searches for gamma-rays from known dwarf galax¬ 
ies [TBUTS] have yielded stringent constraints on the dark 
matter parameter space. They have not yet, however, 
been sufficiently sensitive to cover the full range of cross 
sections favored to explain the Galactic Center excess. 
It has been anticipated that ongoing and planned optical 
surveys will discover a significant number of presently un¬ 
known Milky Way dwarf spheroidal galaxies [T9]-[22]. If 
one of more of these objects happens to be nearby and/or 
contain a high density of dark matter, it could constitute 
an important target for the Fermi Gamma-Ray Space 
Telescope, significantly strengthening their sensitivity to 
annihilating dark matter. 

Very recently, optical imaging data from the Dark En¬ 
ergy Survey (DES) was used to discover nine new dwarf 


galaxy candidates [24]. Shortly thereafter, yet an¬ 
other dwarf candidate (Triangulum II) was discovered 
from within the data from the Panoramic Survey Tele¬ 
scope and Rapid Response System (Pan-STARRS) |25| . 
Of particular interest is the object Reticulum II (also 
known as DES J0335.6-5403) whose proximity (^30-32 
kpc) and spatial extent (half-light radius of 55 ± 5 pc) 
make it likely to be a dwarf galaxy (rather than a glob¬ 
ular cluster) and a very promising target for gamma-ray 
searches for annihilating dark matter. Although spec¬ 
troscopic follow-up will be required to measure the dark 
matter distributions of these systems, it is plausible that 
Recticulum II (or perhaps Triangulum II) could provide 
a gamma-ray signal from annihilating dark matter that 
is brighter than that from any other known dwarf galaxy. 

Two independent groups have already reported the re¬ 
sults of their analyses of Fermi Gamma-Ray Space Tele¬ 
scope data from the directions of Reticulium II and DBS’s 
other new dwarf galaxy candidates. The first of these 
analyses, presented jointly by the Fermi and DES col¬ 
laborations, identified a modest gamma-ray excess from 
the direction of Reticulium II, with a test statistic (TS) 
of 6.7 [26|.^ Although the text of that paper quotes 
a p-value for this excess of 0.06, this includes a trials 
factor of approximately 4, intended to account for the 
range of possible dark matter masses and annihilation 
channels that were scanned over. If one instead employs 
dwarf galaxies to test the dark matter interpretation of 
the Galactic Center excess (for which the spectrum has 
been previously measured), no such trials factor is re- 


^ TS is defined as twice the difference in the global log-likelihood 
between the null and alternative hypotheses. 
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quired. The local significance of this signal (correspond¬ 
ing to p ^ 0.015) is approximately 2.4cr. The other anal¬ 
ysis, carried out by Geringer-Sameth et a/., also finds 
an excess from Reticulium II, and quotes significances 
ranging from approximately 2.3cr to 4.7cr, depending on 
the background model procedure and trials factor that 
are adopted. The local significance (without trials for 
different dark matter models) determined using an em¬ 
pirical sample of nearby spatial regions (as opposed to 
assuming that the background is Poisson distributed) is 
found to be approximately 2.8cr m, or approximately 
0.4cr greater that that found by the Fermi and DES col¬ 
laborations. The excesses identified in each of these pa¬ 
pers is most prevalent at energies between ^2-10 GeV, 
compatible with the spectrum observed from the region 
surrounding the Galactic Genter Gonsidering the 

fact that these two groups apply quite different analy¬ 
sis techniques to this problem, and employ different data 
sets (Ref. [26] utilizes Pass 8 data while the analysis of 
Ref. [27| is restricted to the publicly available Pass 7), it 
is remarkable that their results are so similar. 

In this article, we present our own analysis of the (Pass 
7) Fermi data from the directions of Reticulum II and 
other Milky Way dwarf galaxies with the goal of assess¬ 
ing the characteristics and statistical significance of any 
excess that might exist. We confirm the existence of a 
gamma-ray signal from Reticulum II, and assess the (lo¬ 
cal) significance of this excess to be 3.2cr. And although 
spectroscopic follow-up of Reticulum II will be required 
before this observation can be used to constrain or infer 
the value of the dark matter annihilation cross section, 
for a plausible range of dark matter profiles, this result 
appears to be consistent with dark matter interpreta¬ 
tions of the Galactic Genter signal and the null results 
from other dwarf galaxies. 


II. FERMI DATA ANALYSIS 


In order to calculate the significance of any gamma-ray 
emission observed from a given dwarf galaxy (or dwarf 
galaxy candidate), we examine approximately 6.5 years 
of Fermi-LAT data,^ using the P7REP photons in the 
energy range of 500 MeV to 500 GeV. We exclude events 
arriving with a zenith angle greater than 100°, as well 
as those which do not pass the “Source” photon data se¬ 
lection. We also exclude events that were observed while 
the instrument was not in science survey mode, when the 
instrumental rocking angle was >52°, or when the instru¬ 
ment was passing through the South Atlantic Anomaly. 
For each source, we examine the photons observed within 
a I0°xI0° box centered around the location of the source, 
and divide the photons into lOOxIOO angular bins and 24 
evenly spaced logarithmic energy bins. We analyze the 


2 MET range: 239557417 - 447078115 



FIG. 1: A comparison of the values of the test statistic 

found in our analysis {TShl) to those given in Fermi’s 3FGL 
catalog (TSsfgl)- Overall the agreement between these two 
determinations is good, validating our analysis procedure. 


instrumental exposure throughout the region of interest 
using the P7REP_S0URCE_V15 instrument response func¬ 
tions. We employ the latest model for diffuse galactic 
gamma-ray emission (gll_iem_v05_revl. f it) and the 
latest isotropic emission template for the Source pho¬ 
ton events (iso_source_v05.txt), and include all point 
sources given in the 3FGL Gatalog |28| . 

In our analysis, we follow a prescription as similar 
as possible to that employed by the Fermi-LAT col¬ 
laboration [m [iHi EH- Specifically, we first set the 
global normalization of background sources and the dwarf 
spheroidal over the entire 500 MeV - 500 GeV energy 
range, utilizing the Fermi-LAT gtlike code and the 
MINUIT algorithm. In this phase, we seed the dwarf 
spheroidal spectrum as a simple power-law with an in¬ 
dex of -2.0. We then fix the normalization of all back¬ 
ground components, and employ the pyLikelihood pack¬ 
age to scan the flux of the dwarf spheroidal in each energy 
bin, calculating the delta-log-likelihood {ALG{L)) as a 
function of the source flux. In order to calculate the TS 
for a given dark matter model, we minimize the total log- 
likelihood summed over all energy bins after constraining 
the photon flux by the spectral shape of the dark matter 
model. 

To validate the results of this method, we perform two 
tests. First, we randomly select 100 Fermi-LAT 3FGL 
point sources [28] with |6| > 30°, a “curve significance” 
smaller than 2 (indicating consistency with a power-law 
spectrum) and a TS smaller than 100 in the energy range 
of 300 MeV to 100 GeV. In order to compare our results 
to those given in the 3FGL catalog, we employ the above 
technique with the following modifications. We restrict 
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FIG. 2: The fraction of “blank sky” locations with a test 
statistic (TS) larger than a given value, as empirically deter¬ 
mined for a collection of 1905 randomly selected sky locations 
constrained to lie at a galactic latitude |6| > 30° and at least 
1° (5°) from point-like (extended) 3FGL sources [28|. For the 
blue curve, no additional requirements are placed on the blank 
sky locations. For the red curve, the blank sky locations used 
are additionally required to lie no closer than 0.5° from any 
source listed in the BZGAT, GRATES, GGRaBS, or ATNF 
catalogs (see Sec. IV). The shaded region surrounding each 
curve represents the poisson errors on this determination. In 
generating this figure, we have adopted a spectral shape cor¬ 
responding to a 49 GeV dark matter particle annihilating to 
bb (corresponding to the best-fit mass for the Galactic Genter 
gamma-ray excess my 


our analysis to four years of Fermi-LAT data,^ evaluate 
an energy range of 300 MeV - 100 GeV in 20 energy bins 
utilizing 8 energy bins per decade except for the final bin 
(which was extended to an energy of 100 GeV), and we 
scan the likelihood fits using power-law, rather than dark 
matter motivated, spectral shapes. In Fig. we show the 
distribution of the TS calculated in our analysis (TShl) 
compared to that obtained by the Fermi-LAT collabora¬ 
tion (TS^fgl) in the same energy range. We find that our 
TS values are, on average, slightly (13.5%) lower those 
reported in the 3FGL. We attribute this primarily to the 
fact that we normalize the background by fitting over a 
10° X 10° region, rather than over the entire sky. The 
dashed curve in Fig. [^represents the best-fit gaussian of 
this distribution, with a mean of -0.135 and a standard 
deviation of 0.176. 

Secondly, we apply the “blank-sky” null-test employed 
in previous dwarf spheroidal studies. Specifically, we se¬ 
lect 1905 sky locations with \b\ > 30°, which are 1° re¬ 
moved from any 3FGL source and 5° removed from any 
extended 3FGL source. In this case, we employ the full 
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EIG. 3: The log-likelihood fit of Reticulum II in 24 energy 
bins spanning 500 MeV to 500 GeV. The upper limits corre¬ 
spond to 2cr confidence in each energy range. The white line 
corresponds to the best fit from a 49 GeV dark matter particle 
annihilating to bb. 


6.5 years of data, adopt the default energy range, and 
test the comparison to a 49 GeV dark matter model an¬ 
nihilating to bb (corresponding to the best-fit value of 
the mass for the spectrum of the Galactic Center ex¬ 
cess [9]). In Fig. [^ we show the resulting distribution of 
our blank-sky test locations. While the existence of sys¬ 
tematic errors in the modeling of the gamma-ray back¬ 
ground drives this distribution far from that expected 
from Poisson variations, the result is in good agreement 
with all previous studies. In this figure, we show results 
corresponding to the case in which no additional require¬ 
ments are placed on the blank sky locations (blue), and 
to when the blank sky locations used are further required 
to lie no closer than 0.5° from any source listed in the BZ- 
CAT, CRATES, CGRaBS, or ATNF catalogs (red). This 
will be discussed in more detail in Sec. llYl 


III. RESULTS 

In Fig. [^ we show the delta-log-likelihood {ALG{L)) 
distribution for our analysis of Fermi data from the direc¬ 
tion of Reticulum IF As in both Ref. [26] and Ref. [27], 
we find an excess of events in the bins covering approx¬ 
imately ^2-10 GeV. For a spectral shape corresponding 
to a 49 GeV dark matter particle annihilating to bb (the 
best-fit mass for the Galactic Center excess 0), we find 
a value of TS=17.4 from Reticulum II, corresponding to 
a significance of 3.2cr (see Fig. |^. If we do not impose 
this choice of the dark matter mass, but rather allow the 
mass to float as a free parameter, the value of the TS 
increases only slightly (to 18.1), illustrating the compat¬ 
ibility between this signal and that observed from the 
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logio('^) 

TS (Point-Like) 

TS (NFW-Like) 

Dwarf Name 

Distance (kpc) 

Latitude ('^) 

Ref. [23 

Ref. [30] 

— 49 GeV 

any tudm 

T^DM — 49 GeV 

any tudm 

Reticulum II 

32 (30) 

-49.7 


- 


- 

17.4 

18.1 

- 

- 

Tucana II 

58 (69) 

-52.4 


- 


- 

1.44 

1.82 

- 

- 

Indus I 

69 (100) 

-42.1 


- 


- 

0.0 

0.0 

- 

- 

Horologium I 

87 (79) 

-54.7 


- 


- 

0.09 

0.17 

- 

- 

Phoenix II 

95 (83) 

-59.7 


- 


- 

0.0 

0.55 

- 

- 

Eridanus III 

95 (87) 

-59.6 


- 


- 

0.0 

0.53 

- 

- 

Pictoris I 

126 (114) 

-40.6 


- 


- 

0.0 

0.0 

- 

- 

Grus 1 

(120) 

-58.8 


- 


- 

0.0 

0.40 

- 

- 

Eridanus II 

330 (380) 

-51.6 


- 


- 

0.0 

0.61 

- 

- 

Triangulum II 

30 

-23.8 

- 

- 

0.0 

0.0 

- 

- 

Canis Major 

7 

-8.0 


- 


- 

0.0 

0.0 

- 

- 

Segue 1 

23 

50.4 

19.5 

± 

0.29 

19 36“^°'^^ 
J^y-oo_0.35 

1.07 

1.18 

1.10 

1.72 

Sagittarius 

26 

-14.2 


- 


- 

2.13 

4.33 

- 

- 

Ursa Major II 

32 

37.4 

19.3 

± 

0.28 

iy.nz_o 42 

0.0 

0.32 

0.02 

0.88 

Segue 2 

35 

-38.1 


- 



0.49 

2.09 

- 

- 

Willman 1 

38 

56.8 

19.1 

± 

0.31 

- 

3.94 

4.47 

5.70 

5.89 

Bootes II 

42 

68.9 


- 


- 

0.0 

0.25 

- 

- 

Coma Berenices 

44 

83.6 

19.0 

± 

0.25 

19.02l°;|J 

0.0 

0.0 

0.0 

0.04 

Bootes III 

47 

75.4 


- 


- 

0.0 

0.31 

- 

- 

Bootes I 

66 

69.6 

18.8 

± 

0.22 

18.24l°;« 

0.0 

0.75 

0.09 

0.50 

Draco 

76 

34.7 

18.8 

± 

0.16 

18.84l0;12 

0.04 

0.05 

0.01 

0.01 

Ursa Minor 

76 

44.8 

18.8 

± 

0.19 

-LO-y'J_0.19 

0.0 

1.36 

0.0 

0.99 

Sculptor 

86 

-83.2 

18.6 

± 

0.18 

18 54+0-06 

0.01 

0.02 

1.22 

1.28 

Sextans 

86 

42.3 

18.4 

± 

0.27 

1 7 52+0-28 
i/.OZ_o 18 

0.0 

0.26 

0.0 

0.24 

Ursa Major I 

97 

54.4 

18.3 

± 

0.24 

17 87+0-66 

^ ‘ -0.33 

0.0 

0.07 

0.0 

0.24 

Carina 

105 

-22.2 

18.1 

± 

0.23 

' -0.09 

0.0 

0.03 

0.0 

0.20 

Hercules 

132 

36.9 

18.1 

± 

0.25 

iO.OO_Q 68 

3.09 

4.05 

3.58 

4.11 

Fornax 

147 

-65.7 

18.2 

± 

0.21 

17.83l°;J2 

0.36 

0.94 

1.01 

1.40 

Leo IV 

154 

56.5 

17.9 

± 

0.28 

16.32ll;°e 

0.0 

0.0 

0.0 

0.0 

Canes Venatici H 

160 

82.7 

17.9 

± 

0.25 

17.65l°;« 

0.27 

1.56 

0.50 

1.62 

Leo V 

178 

58.5 


- 


16.37l°J4 

0.0 

0.39 

- 

- 

Pisces H 

182 

-47.1 


- 


- 

0.0 

0.0 

- 

- 

Canes Venatici I 

218 

79.8 

17.7 

± 

0.26 

17 43+0-37 

i / •^o_Q 28 

0.39 

0.47 

0.35 

0.42 

Leo H 

233 

67.2 

17.6 

± 

0.18 

17.97l°;20 

0.0 

0.0 

0.0 

0.0 

Leo I 

254 

49.1 

17.7 

± 

0.18 

17.84l°;20 

0.0 

1.67 

0.0 

1.91 


TABLE I: The distance, galactic latitude, J-factors, and test statistic (TS) of any gamma-ray excess from DES’ nine newly 
discovered dwarf galaxy candidates (top), the new Pan-STARRS dwarf candidate (middle), and the 25 previously known Milky 
Way dwarf galaxies (bottom). Eor each of the DES dwarf candidates, we list distances as reported in Ref. [23] (and as reported 
in Ref. [H]). Eor Triangulum II, we list the distance as reported in Ref. [25]. Eor the known dwarf galaxies, distances are as 
given in Refs. [29l|30]. Although we list only central values for distances, the error bars on these quantities are typically on 
the order of ±(10 — 15)%. The J-factors are averaged over a 0.5° radius around each dwarf, as reported by Refs. [29] and [30] . 
respectively, and are given in units of GeV^cm“^. The TS values listed for each dwarf assume either a spectral shape that 
corresponds to a 49 GeV dark matter particle annihilating to bb (corresponding to the best-fit mass for the Galactic Genter 
gamma-ray excess El) or allowing the dark matter mass to be a free parameter. Eor those dwarf galaxies constrained by stellar 
kinematics, we also show the TS values found when the source is treated as a spatially extended object, corresponding to an 
NEW halo with a scale radius equal to central value reported in Refs. [IHIESI- The most significant detection is from the newly 
discovered and nearby dwarf galaxy candidate Reticulum II. 


Galactic Center. 

In Table [T| we list the TS values found in our analy¬ 
sis for each of the previously known Milky Way dwarf 
spheroidal galaxies, and for the ten newly discovered 
dwarf galaxy candidates. Values are given assuming ei¬ 
ther a spectrum corresponding to the best-fit mass for 
the Galactic Center excess, or for any dark matter mass. 
For a Galactic Center-like spectrum, Reticulum II yields 
the highest significance (TS=17.4), followed by Willman 
1 (3.94), Hercules (3.09), Sagittarius (2.13), Tucana II 
(1.44), and Segue 1 (1.07). Other than that from Reticu¬ 


lum II, no statistically significant excesses are observed. 

Also given in Table |T] are the values of the J-factors 
for each dwarf galaxy with sufficient kinematic informa¬ 
tion (from spectroscopic data) to obtain a determination. 
This quantity is defined as follows: 

j= f [ p^di 

J AO _ J los 

where Aft is taken to be a circle of 0.5° radius around the 
given dwarf, p describes the dark matter density profile 
of the dwarf, and the second integral is performed over 


( 1 ) 
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FIG. 4: The test statistic (TS) as a function of the assumed scale radius of the (NFW) dark matter halo (in degrees, 
tan^s = rs/D) for five selected dwarf galaxies. The blue (green) curves correspond to a spectral shape for the best-fit dark 
matter mass (for ttidm = 49 GeV) annihilating to bb. 


the observed line-of-sight (los). We provide the J-factor 
values as reported by two groups: Martinez et al [29] and 
Geringer-Sameth et al (SD). In general, the dark matter 
profiles of the classical dwarfs are well constrained by stel¬ 
lar kinematics, resulting in relatively well determined J- 
factors. In contrast, the ultra-faint dwarfs (Segue 1, Ursa 
Major II, Willman 1) contain far fewer stars, and exhibit 
much larger J-factor error bars. Deeper measurements, 
capable of detecting more numerous faint stars, will ulti¬ 
mately improve this situation. Although no spectroscopic 
information exists for any of the ten new dwarf galaxy 
candidates, we expect such follow-up measurements to 
occur in the near future. 

For those dwarfs with profiles constrained by stellar 
kinematics, we also list in Table [T| the values of the TS 
found when the source is treated as a spatially extended 
object, rather than as a point-like source. In particular, 
we adopt an NFW-profile for these systems, with a scale 
radius equal to the central value reported in Refs. [ia[29|. 
In Fig. we plot the TS as a function of the halo’s scale 
radius (in degrees, tan6>s = Vq/D) for five of the dwarf 
galaxies under consideration. No strong evidence for (or 
against) spatial extension is observed. The significances 


of Willman 1 and Ursa Major II marginally increase if 
an extended halo is assumed, while the significances of 
Reticulum II and Sagittarius marginally decrease. 

IV. CONTROLLING BACKGROUNDS WITH 

MULTI-WAVELENGTH SOURCE CATALOGS 

In Ref. [31], it was pointed out that Fermi’s sensitivity 
to dark matter annihilation in dwarf spheroidal galaxies 
could be increased by taking into account information 
available in multi-wavelength source catalogs. In partic¬ 
ular, a significant fraction of the highest TS points in 
the “blank sky” correspond to the locations of unresolved 
blazars, radio galaxies, and star forming galaxies. By 
making use of only regions of the “blank sky” which are 
not near sources listed in multi-wavelength catalogs, it is 
possible to reduce the contamination from such sources. 

In Fig. the TS distribution of the high-latitude 
blank-sky is shown without utilizing multi-wavelength in¬ 
formation (blue), and after avoiding all locations located 
within 0.5° of any source listed in the Roma-BZCAT 
Multi-Frequency Catalog of Blazars (BZCAT) [32], the 
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Combined Radio All-Sky Targeted Eight-GHz Survey 
(CRATES) catalog [33], the Candidate Gamma-Ray 
Blazar Survey (CGRaBS) catalog [31|, or the Australia 
Telescope National Eacility (ATNE) pulsar catalog [35| 
(red). The application of this cut significantly reduces 
the fraction of the sky with large TS values. 

To take this multi-wavelength information into ac¬ 
count, we apply the following procedure in our analysis. 
Eor a given dwarf galaxy (or dwarf galaxy candidate), we 
check the catalogs described in the previous paragraph 
for any sources located within 0.5°. If any are found, we 
re-run our analysis, including in the background model 
a source at that location. We then take the new TS of 
the dwarf, and see how many locations on the blank sky 
yield a higher TS, when a background source is included 
at the location of the nearby catalog source. We then use 
the fraction of high-TS blank-sky locations to calculate 
the p-value and significance of any dwarf galaxy excess. 

This procedure is most important in the case of Retic¬ 
ulum II, which is located 0.44° from the source CRATES 
4033553-543025.^ Given the large number of sources con¬ 
tained in these catalogs, this is not a particularly sur¬ 
prising (we estimate a probability of ^20% that at least 
one source would reside within 0.5°). The analysis with 
the extra background source at the CRATES location 
resulted in a TS of 13.4 from Reticulum II (a modest 
reduction from 17.4).^ After re-running our analysis on 
all “blank-sky” locations with TS>13.4, including in the 
background model sources at the locations of the multi¬ 
wavelength catalog sources, we find that only 3 out of 
1905 of the blank-sky locations yielded a more significant 
excess. This corresponds to a p—value of 0.001575 and a 
detection significance of 3.2cr (compared to 3.0 <j, which 
is found if no multi-wavelength information is utilized). 
If we include the poisson error bars around the number 
of 3 blank-sky locations, the corresponding significance 
covers the range of 3.0 to 3.4cr. 

As we were finalizing this paper, it was pointed out 
to us that the faint radio source PMN 40335-5406 in the 
Parkes-MIT-NRAO catalog is located 0.105° from the 
location of Reticulum II.^ The PMN catalog contains a 
much larger number of sources (50,814) than the catalogs 
employed in our analysis, and it is not particularly sur¬ 
prising that a such source resides near Reticulum II (or 
near any other location in the high-latitude sky). Specifi¬ 
cally, at the declination of Reticulum II, the PMN survey 
contains approximately 17 sources per square degree, cor¬ 
responding to a 17% chance that a catalog source exists 


^ The presence of this source was pointed out to us by Eric Carlson. 

^ We note that the spectral shape absorbed by the CRATES source 
is very hard, and quite unlike that of the gamma-ray emission 
from typical radio sources. We consider it unlikely that this 
source contributes significantly to the gamma-ray fiux observed 
from the direction of Reticulum II. 

® We thank Alex Drlica-Wagner for bringing this source to our 
attention. 



Right Ascension (J2000) 

FIG. 5: The test statistic (TS) over the region of the sky 
surrounding Reticulum II. Although the excess gamma-ray 
emission is localized to a region of approximately ~ 0.1° ra¬ 
dius, it is not possible to spatially discriminate Reticulum 
II (black square) from the location of the radio source PMN 
40335-5406 (blue oval). 

within a 0.1° ROI around Reticulum II. In Eig. we 
plot the value of the TS found by our analysis for re¬ 
gions of the sky in the vicinity of Reticulum II. Although 
we can localize the excess gamma-ray emission to a re¬ 
gion of approximately ^0.1° radius, it is not possible to 
spatially discriminate Reticulum II (black square) from 
the location of the PMN source (blue oval). The PMN 
source exhibits continuous radio emission that is consis¬ 
tent with a blazar origin. When re-running our analysis 
for a power-law gamma-ray spectrum (as is observed from 
blazars) we find a best-fit spectral index of 2.04, yielding 
TS=13.1. In comparing this to the TS value found earlier 
in this study, we conclude that a dark matter-like spec¬ 
trum is preferred to a blazar-like power-law at the level of 
approximately ^ 2cr. Spectroscopic observations of this 
source would be valuable, as they could aid in determin¬ 
ing its nature and inform us as to its likely gamma-ray 
luminosity and spectrum. 

V. A SELF-CONSISTENT INTERPRETATION 

In this section, we consider the excess observed from 
Reticulum II, along with the lack of significant detections 
from other Milky Way dwarf galaxies, and ask whether 
these results are mutually consistent. Focusing on the 
case in which annihilating dark matter is responsible for 
the Galactic Center excess (moM = 49 GeV, for the 
case of annihilatiions to 66), the lack of significant ex¬ 
cess emission from the known dwarf galaxies constrains 
cru < 1.3 X 10“^^ cm^/s [16]. And while this constraint 
is compatible with dark matter interpretations of the 
Galactic Center excess, the normalization of the Galac- 
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tic Center signal implies that the cross section is unlikely 
to be smaller than this value by more than a factor of a 
few. More specifically, if we allow the overall normaliza¬ 
tion and the scale radius of the Milky Way’s dark matter 
halo profile to vary within the range allowed by dynami¬ 
cal constraints (piocai = 0.24 — 0.46 GeV/cm^, = 8 — 35 
kpc |36l [37]), we find consistency with an annihilation 
cross section as small as ^ 3 x 10“^^ cm^/s. From this 
perspective, the prospects for the future detection of a 
gamma-ray signal from one or more dwarf galaxies ap¬ 
pears encouraging. 

For an annihilation cross section at the upper limit of 
Ref. [16] {(TV ^ 1.3 X 10“^^ cm^/s), the normalization 
of the signal from Reticulium II requires approximately 
log J ^ 19.6-20.1. Noting the empirical (and approx¬ 
imate) relationship between the distances and J-factors 
of ultra-faint dwarfs. Ref. [26] points out that Reticulium 
II might be expected to have a somewhat smaller value, 
log J ^ 19.3, although even a value as high as of 20.1 
would not be a particuarly significant outlier. The ne¬ 
cessity of a large J-factor for Reticulum II (if its gamma- 
ray excess is from annihilating dark matter) can also be 
seen from the results of our analysis, as shown in Ta¬ 
ble [Ij Roughly speaking, the predicted value for the TS of 
a given dwarf is proportional to its gamma-ray flux, and 
thus to its J-factor. The modest TS values observed from 
Segue 1 and Ursa Major II suggest significantly lower J- 
factors for these systems than for Reticulum II. Given 
this situation, we eagerly await the spectroscopic follow¬ 
up of Reticulium II. If the gamma-ray excess from this 
source in fact originates from annihilating dark matter, 
we should anticipate a large value for its J-factor, likely 
in excess of GeV^ cm“^. 


VI. ARE STATISTICAL FLUCTUATIONS 
SUFFICIENT TO EXPLAIN THE DIFFERENCES 
BETWEEN THE RESULTS FOUND USING PASS 
7 AND PASS 8 DATA? 

At face value, our determination of TS=I8.I from the 
direction of Reticulum II appears to be in conflict with 
the more modest value of 6.7 quoted by the Fermi Col¬ 
laboration [26]. The most significant difference between 
these two analyses is in the data sets that are being con¬ 
sidered: our analysis makes use of the publicly avail¬ 
able Pass 7 data, whereas the Fermi Collaboration pa¬ 
per utilizes the more recent Pass 8 data set. Intrigu- 
ingly, a similar discrepancy can be seen in a compar¬ 
ison of the Fermi’s Collaboration’s Pass 7 and Pass 8 
studies of known dwarf galaxies. In particular, whereas 
the Fermi Collaboration’s Pass 7 analysis revealed an ex¬ 
cess from three ultra-faint dwarfs (Segue I, Ursa Major 
II, and Willman I), at a level of TS ^ 10 [18], no ex¬ 
cess was observed according to their more recent Pass 
8 paper m- This would-be conflict between the Fermi 
Pass 7 and Pass 8 dwarf papers is somewhat surprising in 
light of the fact that these two analyses make use of data 


taken over significantly overlapping time periods (over 
the first four years of Fermi’s mission; the more recent 
analysis adds two more years of data to this set). The 
Fermi Collaboration points out, however, that after tak¬ 
ing into account the new event selection associated with 
the transition from Pass 7 to Pass 8, the overlap between 
these two data sets is not particularly large; only ^30% 
of the 1-10 GeV photons used in the most recent analysis 
were also employed in their earlier study (put another 
way, approximately 65% of those events in Fermi’s Pass 
7 dwarf analysis were also included in the recent Pass 8 
study) [16]. Assuming that the previous TS ~ 10 excess 
was the result of a statistical fluctuation, we estimate 
that there was an approximately ^10% chance that the 
more recent data set would yield TS < 1 (after taking 
into account differences in effective area and exposure 
time). So in this respect, we concur with the conclusion 
of the Fermi Collaboration. 

In order to compare our results directly to those 
of Ref. [18], we reanalyzed the three ultra-faint dwarf 
spheroidal galaxies (Segue 1, Willman 1, and Ursa Ma¬ 
jor II) using the P7V6 Reprocessed data, employing 
in this test only four years^ of data and utilizing only 
those events which pass the P7V6REP “Clean” selection 
cut (and the gll_iem_v05 .f it and iso_clean_v05.txt 
models for the diffuse and isotropic emission, respec¬ 
tively), and examine the data a in 14° x 14° region 
around each ultra-faint dwarf in 0.1° spatial bins. As¬ 
suming a spectrum corresponding to a 25 GeV dark mat¬ 
ter particle annihilating to 56, we find that these three 
dwarfs acquire TS values of 1.18, 2.43, and 2.43, respec¬ 
tively. If we model these sources as extended objects 
(NFW-like) rather than as point sources, we find TS val¬ 
ues of 0.89, 3.05, and 2.73. In either case, our TS values 
are similar to but somewhat lower than the TS ~ 10 re¬ 
ported in Ref. m- Comparing this to the results shown 
in Table |T] however, we see little indication that the addi¬ 
tion of 2.5 years of data has reduced the overall TS from 
these three ultra-faint dwarfs (Ursa Major II’s TS fell, 
while that of Willman 1 increased by a similar amount). 
In light of this somewhat confusing situation, we eagerly 
await the public availability of Pass 8 data. 


VH. SUMMARY AND CONCLUSIONS 

In this article, we have revisited the gamma-ray emis¬ 
sion from known Milky Way dwarf galaxies, and from 
the dwarf galaxy candidates recently discovered in the 
data from DES and Pan-STARRS. Of particular interest 
are the new dwarf candidates Reticulum II and Triangu¬ 
lum II, which are each located at a distance of only ^30 
kpc from Earth, making them promising targets for dark 
matter searches. Our analysis of Eermi data from the di- 
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rection of Reticulum II identifies an excess of gamma-rays 
with a local statistical significance of 3.2cr. This is slightly 
higher than, but not dissimilar to, that reported in the 
previous studies of other groups [Ml [27] . We also confirm 
that Reticulum IPs gamma-ray excess is most prominent 
at energies between ^2-10 GeV, in good agreement with 
the spectral shape of the excess previously reported from 
the region surrounding the Galactic Genter. We do not 
observe any significant 7 -ray emission from the direction 
of Triangulum 11. 

Looking forward, spectroscopic follow-up of Reticulum 
II and the other new dwarf galaxy candidates will be im¬ 
portant for interpreting this data. In order for this excess 
to be compatible with the lack of significant gamma-ray 
detections from other dwarf galaxies (most importantly. 
Segue 1 and Ursa Major II), Reticulum II must contain a 
high density of dark matter, corresponding to J > 
GeV^crn”^. A measurement of Reticulum IPs J-factor 
that is much smaller than this value would place seri¬ 
ous doubt as to any dark matter interpretation of its 


excess. Additional data from Fermi will also have much 
to bear on this question. With 50% more data, such as 
could be acquired over the next few years, we estimate 
that the detection of Reticulum IPs gamma-ray emission 
could exceed TS=25, corresponding to approximately 4a 
significance, and comparable to the threshold for mem¬ 
bership in Fermi’s point source catalogs. 
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